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Abstract. Relativity theory is useful for understanding the phenomenology of the magnetoelectric effect 
of the antiferromagnet chromium sesquioxide Cr203 in two respects: (i) One gets a clear idea about the 
physical dimensions of the electromagnetic quantities involved, in particular about the dimensions of the 
magnetoelectric moduli that we suggest to tabulate in future as dimensionless relative quantities; (ii) 
one can recognize and extract a temperature dependent, 4-dimensional pseudoscalar from the data of 
magnetoelectric experiments with Cr2 03. This pseudoscalar piece of Cr2 03 is odd under time reflections 
and parity transformations and is structurally related ("isomorphic") to the gyrator of electric network 
theory, the axion of particle physics, and the perfect electromagnetic conductor of electrical engineering. 

PACS. 75.50.Ee Antiferromagnetics - 03.50.De Classical electromagnetism. Maxwell equations - 06.20.F- 
Units and standards - 46.05.+b General theory of continuum mechanics of solids - 14.80.Mz Axions and 
other Nambu-Goldstone bosons 



1 Introduction 

Magnetism is caused by moving electric charge, and the 
motion of an object is closely related to relativity theory; 
at least if the charge moves with a speed v that cannot 
be neglected as compared to the speed of light c. Accord- 
ingly, in the theory of magnetoelectricity that deals with 
materials where a magnetisation M is induced by an elec- 
tric field E 01 a, polarization P by a magnetic field B (or 
a magnetic excitation H), magnetic and electric fields are 
closely intertwined, and it is wise to put the electrodynam- 
ical theory in a relativistic form right away; and this all 
, the more since Maxwell's theory is known to be a relativis- 
tic theory covariant under all coordinate transformations 
(see Post [46]). 



2 Physical dimensions in magnetoelectrics 
2.1 Permittivity and permeability 

Before one formulates a physical theory, in the case un- 
der consideration electrodynamics of material media, one 
has to make a basic decision: Do we want to formulate 
the laws in terms of "quantity equations" or in terms of 
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"numerical equations"? A quantity A ~ {A}[A] is always 
a "numerical value" {A} times a "dimension" [A]; and 
quantity equations interrelate these quantities. By their 
very definition quantity equations are valid in any system 
of units. Thus, if one puts the Maxwell equations in such 
a form, they are valid in SI (the international system of 
units), in the Gauss system of units, etc. By contrast, in 
the book of Jackson [29] , for instance, on the top of each 
second page, it is indicated for which system of units (SI 
or Gaussian) the equations displayed are valid. Clearly, 
this distracts from the universal nature of the Maxwell 
equations — and we will not follow that convention. Also 
the elementary particle physicist prefer to formulate the 
laws in such a way that they are only valid in a certain 
system of units, with h — 1, c = 1, etc. 

In order to recognize the archetypal structure of the 
Maxwell equations, one should write them, as already Max- 
well has done himself — who also discusses physical dimen- 
sions in quite some detail — in terms of quantity equations, 
see Post [46] or [18]. We will follow Post in this context. 

In the past these matters of dimensions and units caused 
heated debates. One can get a good idea of this from the 
article of Brown [9] . He gives a good and complete descrip- 
tion of how to change from one system of units to another 
one. However, in contrast to Post, he argues that choice of 
a dimension is more or less a convention. We disagree with 
this position. In our opinion the dimension of a quantity 
encodes the operational procedure of how to measure this 
quantity. 
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The dimension of a quantity is an intrinsic and defin- 
ing property of a quantity. We cannot, by suitable units, 
kill the dimension of a physical quantity. A dimension 
is something invariant. The electric excitation D, for in- 
stance, can be measured in vacuum with the help of the 
Maxwellian double plates implying the dimension [D] — 
charge/area = q/i"^, where q denotes the dimension of 
charge and £ the dimension of length. Similarly the electric 
field strength [E] = force/charge = {w/£)/q = {U It / (.) / [It) 
= U /£, with U, I, and t as dimensions of work (energy), 
voltage, current, and time, respectively. Clearly then, quite 
independently of the system of units chosen, the permit- 
tivity defined by the ratio of the electric excitation and 
the electric field e = D/E has the intrinsic dimension 

[s] = [D]/[E] = {q/e^)e/u^^ iit/£^)e/u = i/{Rv), 

where R and v are the dimensions of resistance and veloc- 
ity, respectively. In particular, also in vacuum [eq] = Y/v^ 
with y = as the dimension of admittance. Note that 
we talk about dimensional analysis, not about a system of 
units. Our discussion looks like one relating to SI; however, 
this is not the case — rather SI also uses quantity equations 
and hence some equations in SI resemble those we discuss 
here. In the Gaussian system, for example, one introduces 
effectively a new quantity d :— D/sq and thereby finds 
in vacuum d — E. However, d carries no longer the di- 
mension of an excitation and, accordingly, is only some 
superficial notational construct without operational inter- 
pretation. The conventional vacuum relation D — E, used 
in Gaussian units, defies a reasonable dimensional analy- 
sis. 

For ordinary linear and isotropic materials we can write 
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where H is the magnetic excitation and D the electric ex- 
citation that are induced by the externally applied fields E 
and B. We immediately recognize that the moduli charac- 
terizing the material, namely the (relative) permittivity s-^ 
and the (relative) permeability fx^ are both dimensionless, 
quite independently of any system of units. This conve- 
nient parameterization is achieved as soon as we introduce 
some standard and constant Eq and /ioj with the dimen- 
sions of Y/v and l/{Yv), respectively. 

It turns out in a relativistic analysis of the Maxwell 
equations that an admittance — or resistance as its reci- 
procal — are four-dimensional scalar quantities even in the 
Riemannian spacetime of general relativity. This is why we 
concentrate on such a quantity with a dimension of ad- 
mittance. That a velocity is inherent in electrodynamics 
is well-known and was already established by the experi- 
ment of Weber-Kohlrausch before the advent of Maxwell's 
theory, see [S]. This vacuum speed of light c is a special- 
relativistic 4-dimensional scalar as was uncovered when 
Minkowski set up a manifestly Poincare-covariant version 
of electrodynamics (compare JJIj). Therefore, two 4-dimen- 
sional scalar quantities, namely the vacuum impedance Yq 
and the vacuum speed of light c characterize the vacuum 



from an electrodynaniical point of view. In different sys- 
tem of units they have different values, in SI, for instance, 
we have Yq = ^Eq/^o ~ 1/(377]?) and c = l/y^eo7*o ~ 
3 X 10* m/s. But the dimensionless permittivity/permea- 
bility, £r and /ir, keep their values in all systems of units. 



2.2 Four-dimensional excitation and field strength 
tensors 

In order to extend this characterization of materials by 
dimensional quantities also to magnetoelectric moduli, we 
have first to turn to the 4-dimensional representation of 
electrodynamics. Each of the two sets {D, H) and {E, B) 
constitutes a relativistic "couple" that is, we can define 
4-dimensional 2nd rank antisymmetric tensors according 
to [\v = 0,1,2,3) 



F - {Fx,) - -{F,x) 



D^' 
-Z?i i/3 -H2 
-D^ -H3 Hi 
-D^ H2 -Hi 

^ ~Ei —E2 — -E3 ' 
El -B'^ 
E2 -B^ B^ 
\ E3 B^ -B^ 



(2) 



(3) 



Even if one is not familiar with these expressions, we 
can see by a dimensional analysis that each couple is 
really intrinsically interrelated. We have [H] — Ijl — 
ql{tl), \D\ = qie, and \D\ = {\lv)\H\. Accordingly, 
D and H are both charge-related quantities distinguished 
only by the dimension of velocity. This dimension of ve- 
locity mediates between time and space components like 
for a coordinate where \x^\ = i, \x^\ = I = vt^ that is, 
[a;°] = (\lv){x\ . . . Analogously, we have [(S"^] = {D^\ = 
{l/v)[Hi] = (l/w)[(5'3] etc. 

For the couple {E, B) we have, with as dimension of 
magnetic flux, [E] = U/i = (t)/{t£) and [B] = 0/^^ that 
is, [E] = v[B]. In other words, E and B are flux-related 
(or, via the Lorentz force, force-related) quantities. 

Apparently, the matrices in ([2]) and ([3]) can be consis- 
tently set up from a dimensional point of view, including 
the factors involving the velocity, provided the mentioned 
couples are used. Thus, [(S] ^ charge, [F] ^ flux, and 
[&] X [F] ^ charge x flux = action; for a detailed discus- 
sion see [IH]. We recognize thereby that the scalar (den- 
sity) with the dimension of an action, 
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(4) 



must have a fundamental importance in electrodynam- 
ics: It represents the Lagrangian of the electromagnetic 
fleld and is equally well a 4-dimensional general relativis- 
tic scalar. 

In the literature these results are hardly recognized. 
Landau-Lifshitz [Mj, for example, like most of the engi- 
neers, see Sihvola and Lindell [SS], take the bastard pair 
{E, H) for their considerations. From a dimensional as 
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well as from a relativistic point of view, there is no legit- 
imacy to give birth to such a hybrid. All dimensions in 
relations get mixed up by such an unholy pairQ 



2.3 Dimensionless magnetoelectric moduli 

According to ([2]) and ([3]) , the general local and linear con- 
stitutive law is expected to be of the general form (S oc F. 
The simplest case is that of an isotropic medium. If then, 
in the magnetoelectric case, a magnetic field B (or an 
electric field E) is applied, the additional response will be 
an electric excitation D (or a magnetic excitation H): 



formulate, following Tamm [63], Post [46], and our dis- 
cussion in [i21 , a 4-dimensional covariant law for a local, 
linear, and anisotropic medium thereby generalizing the 
corresponding isotropic law (l6|),([7]). 
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= r, H = [/3,] = l 
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Since the dimensions of ao and /3o are the same, we 
can put them equally to the vacuum admittance: ao — 
Po = Yq. If we allow an electric and a magnetic field to be 
present at the same time, then ([T]) combined with lU for 
isotropic media leads to the ansatz 



where the constitutive tensor x of rank 4 and weight -1-1 
is of dimension [x] = Y; thus, [^] = 1. Because of the an- 
tisymmetry of (S'^'' and F^-k, we have t^^'"^'* = _t^Ai/k(t _ 
_^i^\aK^ An antisymmetric pair of indices corresponds, 
in four dimensions, to six independent components. Thus, 
the constitutive tensor can be considered as a 6 x 6 matrix 
with 36 independent components. 

A 6 X 6 matrix can be decomposed in its tracefree sym- 
metric part (20 independent components), its antisymmet- 
ric part (15 components), and its trace (1 component). On 
the level of x^'^'^'^j this decomposition is reflected in 
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D = Yq[ er^E + a,B 
H^Yol —cB + f3,E 



Mr 



(6) 
(7) 



[Yq] = admittance , 



[c] = velocity, 

m - 1 . 



Note that Yq/c — Eq and Yqc = l/^o- 

Hence we see that all electric, magnetic, and magne- 
toelectic moduli, independent of any system of units, can 
be described by the dimensionless quantities £i-, fj,,-, a^, Pr- 
This has been shown, amongst others, by Post [46], O'Dell 
[H], Ascher [T|[^. and Hehl et al. [11]. Units such as 
those used in Ref.[16', p. R128, for example, namely (V/cm)/Oe 
should be abandoned. They complicate the comparison of 
different experiments. 

If experimentalists collect magnetoelectric constants, 
then our suggestion is to tabulate them in their dimen- 
sionless relative form, since these values are totally inde- 
pendent of the system of units used and represent their 
genuine values. 



The third part, the axion part, is totally antisymmet- 
ric and as such proportional to the totally antisymmetric 
Levi-Civita symbol, ;= a'e^'"^'^. We remind our- 

selves that the Levi-Civita symbol e'-*""^'* = -|-1 or = — 1 
depending whether Xvan denotes an even or an odd per- 
mutation of the numbers 0123, respectively; it is zero oth- 
erwise, see [61j. We denote pseudoscalars with a tilde. 
The second part, the skewon part, is defined accord- 
(2)^/x,.Ap ._ i(^M^Ap _ ^Appi.^ ^.jjg constitutive 



mg to y^'x'^'^^^' := ^[x " - X 

equation can be derived from a Lagrangian, which is the 
case as long as only reversible processes are considered 
then (2);^A^'^« 



symmetries Wx^"""^ = W^ctkau 



A 

constitutive relation now reads 



0. The principal part W-yA"'^'^ fulfills the 



and (1)^'^ 



= 0. The 
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2.4 Local and linear magnetoelectric constitutive law 

Now that wc have an elementary understanding of the 
dimensions of the magnetoelectric moduli, we proceed to 

^ There is one situation, however, in which the non- 
relativistic pair {E, H) can be useful. Since in the Maxwell 
equations only D and B enter as time derivatives, see curl H = 
j + D and curl E + 13 = 0, the pair {D,B) and by implication 
{E,H) are useful for formulating the initial value problem in 
electrodynamics, as was pointed out by Perlick |45) . In this 
situation one has to leave the manifestly relativistic covariant 
formulation of electrodynamics and introduces necessarily a 
particular slicing of the 4-dimensional spacetime that destroys 
the manifest covariance of the theory. 



with Si- as dimensionless, [a^] — 1. 

In order to compare (jlOp with experiments — the law 
([6]),([7]) is a special case of it — we have to split (|TOl) into 
time and space parts. As shown in j20j in detail, we can 
parameterize the principal part by the 6 permittivities 
^ab _ ^ba^ -j-j^g g permeabilities ^afc — l^bai and the 8 mag- 
netoelectric pieces (its trace vanishes, = 0, sum- 
mation over c!) and the skewon part by the 3 permittivi- 
ties ria, the 3 permeabilities m", and the 9 magnetoelectric 



pieces Sq . Then ([T 



D" = is" 
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can be rewritten as (a, 6, c = 1, 2, 3) 
Eb 

aB", (11) 
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- Sb 



+ (- 
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')Ek-aEa. (12) 



Here e°'"^ = Eabc = ±1, are the 3-dimensional Levi-Civita 
symbols. As can be seen from our derivation, a = Si-lo is a 
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4-dimensional pseudo- (or axial) scalar, whereas Sc'' is only 
a 3-dimensional scalar. The cross-term 7";, is related to the 
Fresnel-Fizeau effects. The skewon contributions m'^,nc 
are responsible for electric and magnetic Faraday effects, 
respectively, whereas the skewon terms Sa'' describe opti- 
cal activity. Equivalent constitutive relations were formu- 
lated, amongst others, by Serdyukov et al. |55j . by Lindell 
[55] , and by Spaldin et al. [62j. Magnetoelectric effects in 
spiral magnets and in ferromagnets were recently studied 
by Mostovoy '39] and Dzyaloshinskii [14], respectively. 

From ([II]),([ni), for = where are the com- 
ponents of the 3-dimensional (contravariant) metric, {/Li~°)ah 
= {^~^)gab, and all other magnetoelectric moduli vanish- 
ing, we recover the special case 

D" = eE" +55", (13) 

Ha - (M"')^a - aEa , (14) 

with = g'^^'Eb and Ba = gabB'' . This is the isotropic 
case of (HI),©; but additionally, we recognize that in ([7]) 
we have to require 

/?r = -«r (15) 

for consistency. 

Incidentally, if necessary, the linear law ([8]) can be 
amended with a quadratic term according to the scheme 

©^"^ = \ X^'""'F,, + ^r,^-""^^^ F^^F^r , (16) 

with an additional constitutive tensor 77 of 126 indepen- 
dent components. In the reversible case they reduce to 56 
independent components. 



To (ii): If a is spacetime dependent, a = a{x), then 
it emerges in the field equations of electrodynamics, as 
had already been shown by Wilczek [70] in his axion- 
electrodynamics. In the constant case, the a jumps at the 
interface between the medium and vacuum and yields a 
contribution to the corresponding boundary conditions, 
see IH]. Moreover, at least since the very exact measure- 
ments of Wiegelmann et al. [SH] no doubt was possible 
at the violation of the Post constraint in the case of the 
antiferromagnet Cr203. 



3.2 Dzyaloshinskii's theory for the magnetoelectric 
effect of the antiferromagnet Cr2 03 



Landau and Lifshitz (1956, see the first Russian edition of 
[34]) predicted the magnetoelectric effect as a phenomenon 
that "results from a linear relation between the magnetic 
and the electric fields in a substance, which would cause, 
for example, a magnetization proportional to the elec- 
tric field..." It "can occur for certain magnetic crystal 
symmetry classes." Dzyaloshinskii [IT], on the basis of 
neutron scattering data and susceptibility measurements, 
found out that the antiferromagnetic chromium sesquiox- 
ide Cr203 has the desired magnetic symmetry class. Start- 
ing from a thermodynamic potential quadratic and bilin- 
ear in E and H, he developed a theory for the magneto- 
electic constitutive relations for Cr203. Written as quan- 
tity equations, they are valid in an arbitrary system of 
units and read 



3 The relativistic pseudoscalar a 

3.1 Violation of the Post constraint 

According to some hand- waving arguments of Post, see 
[46] . p. 129, the pseudoscalar a in ([TT|) . ([T2)) (and even its 
first derivative) should vanish for the vacuum as well as for 
all media. Lakhtakia, see [32] j dubbed the relation a = 
the "Post constraint" and pushed it as a presumed con- 
sequence of Maxwell's equations quite vigorously. Lakhta- 
kia's arguments have never been convincing to a deter- 
mined minority, see, for instance, Sihvola and Tretyakov 
[561160]. More recently Lakhtakia ^33] admitted that he 
was driven to the Post constraint by the following two 
prejudices [33] : 

"(i) The idea of a nonreciprocal but isotropic medium 
appears oxymoronic [contradictory]." 

"(ii) The idea of a constitutive parameter that van- 
ishes from the macroscopic Maxwell equations for a linear 
medium but leaves behind its spatiotemporal derivatives 
appears bizarre." 

To (i): This is the case of ([13]) , ([H]) . If a material is 
"charged" with the 4-dimensional pseudoscBlai a — admit- 
tedly a case that is rare in nature — then an orientation 
dependence of a is implied; why should this contradict 
the 3-dimensional isotropy of the medium? 
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(17) 
(18) 
(19) 

(20) 
(21) 
(22) 



The z-axis is parallel to the optical axis of Cr203. Here 
we have (relative) permittivities parallel and perpendicu- 
lar to the z-axis of the crystal, namely e||,e_L, analogous 
(relative) permeabilities ^J.± and magnetoelectric mod- 
uli a||, a±. 

As we can see from pT]) and (|12l) . we have to get 
(D, H) on the left hand side and (E, B) on the right hand 
side in order to end up with a constitutive law that is 



Friedrich W. Hehl et al.: Magnetoelectric Cr203 and relativity theory 



5 



"written in a relativistically covariant form. We find, 



D^ = Yo 



Dy=Yo 



= Yr 




and 
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(23) 
(24) 

(25) 

(26) 
(27) 
(28) 



We can compare these equations with (ITT|) and p2|) . 
Since Dzyaloshinskii assumed reversibiUty, the skewon piece 
with its 15 independent components vanishes identically. 
Hence (ITlT) and p2|) reduce to 



Ha = ti;jB'' -J^Eb-SEa, 



(29) 
(30) 



with 21 independent moduli. The 4-dimensional pseudo- 
scalar a (alternatively called the axion parameter) repre- 
sents 1 component. We compare ([29]) and ([30]) with ((23] 
to and note that Dzyaloshinskii used Cartesian coor 
dinates such that B'^ = B^ etc. Then we find the permit- 
tivity 



_a6 



>0 

c 



Ml 
ei 
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(31) 



the impermeability 



' 
Yoc\ A*:' 

^i7\ 



(32) 



the magnetoelectric 7 matrix 

3 VMi Mil/ \o 0-2, 



(33) 



We called it relativistic, since it occurs in the context of 
the relativistic (E, B) system. Later, the general form of 
these matrices for the crystal structure of Cr203 was con- 
firmed by O'Dell |31], amongst others. Note that e||, e^, /.t||, 
/.tj_, a||, and a_L, according to their definitions (fT7|) to (|^ . 
are measured in an external E and/or an external H field. 

3.3 Experiments of Astrov, Rado & Folen, and 
Wiegelmann et al. 

The magnetoelectric effect for Cr203 was first found ex- 
perimentally by Astrov [3] for the electrically induced mag- 
netoelectric effect (called MEe in future) and by Rado & 
Folen Wj] for the magnetically induced magnetoelectric ef- 
fect (MEh). In both investigations single crystals of Cr203 
were used. In the MEe experiments [4l[47]. Eqs. (l20l) to 
were verified [H switched off) and in the MEh experi- 
ments 47 Eqs.in]) to IH]) [E switched off). In particular, 
Rado & Folen made both type of experiments and found 
that the magnetoelectric moduli a±_ and a\\ for MEe ex- 
periments coincide with those of the MEjj experiments. 
This proves the vanishing of the skewon part i'^)y^^'^<^'^ of 
the constitutive tensor ■^^'^'^'^ for Cr203. 

Accordingly, these experiments confirmed Dzyaloshin- 
skii's theory for Cr203 below the spin-flop phase. Further 
experiments were then done mainly for the MEh case. Par- 
ticular accurate measurements are those of Wiegelmann et 
al. [55] J see also [57115^ . They took magnetic fields B as 
high as 20 tesla and measured from liquid Helium up to 
room temperature. Wiegelmann et al. took a quasi-static 
magnetic field and thereby proved explicitly that mea- 
surements with magnetic fields of some kilo hertz can be 
extrapolated to static measurements. The sign of a_L(T) 
relative to Qf|| (T) was left open in [68] • Hence we took that 
from Astrov [4]. 

We plotted the values of ^^(r) and a\\{T) from the 
references quoted. We have shown in [21j that the perme- 
abilities « /^ii « 1. Then we can determine a by means 
of (|34l) . The values we found [21] are displayed in Fig. 1. 
Up to about 163 K, the pseudoscalar a is negative, for 
higher temperatures positive until it vanishes at the Neel 
temperature of about 308 K. Its maximal value we find at 
285 K: 



5„,ax (at 285 K) 



3.10 X lO"'' Yo. 



(36) 



Consequently, we have definitely a nonvanishing pseudo- 
scalar a — cirYo with a maximal Sr of the order of 10^'*. 
It is a small effect, but it does exist. 



and the pseudoscalar (or axion) piece 



M-L 



Mil 



Yo. 



(34) 



Conventionally, in the magnetoelectric literature the 
7-matrix and 5 are collected in the "relativistic" matrix 



a b 1 b + adb = Yo 




(35) 



3.4 New predictions 

3.4.1 Search for the first cubic magnetoelectric crystal 

As we saw in Sec. 13.31 the parallel and perpendicular per- 
meabilities for Cr203 are approximately equal, M|| ~ M-L = 
fi. Then the magnetoelectric 7 matrix of (|33p becomes 



7 b 



Yo a_L — «! 



Mr 




(37) 
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100 125 150 175 200 225 
TEMPERATURE T, K 



Fig. 1. The pseudoscalar (or axion) piece a of Eq. (|34|l of 
Cr2 03 in units of the vacuum admittance Yo as a function 
of the temperature T in kelvin; in SI, Yo « 1/(377 ohm); for 
details see |211. 



The question is now: Can we find a substance in which 



a± 



(38) 



that is, in which the matrix j°'b vanishes for all temper- 
atures ? This would be interesting to know since then one 
would have a substance, like in (fT3)) and (fT4)) . in which the 
only magnetoelectric piece would be the pseudoscalar (or 
axion) piece a. 

We can take from Rivera [35] j Fig-2, box 16, that three 
equal entries in the diagonal of the matrix ([37]) are only 
possible in the following five antiferromagnetic cubic point 
groups: 

23, m'3', 432, 4'3m', m'3'm' . (39) 

Schmid [53], p. 23, has explained that the Cr-Br- and Cr-I- 
boracites are antiferromagnetic and cubic below 17 K and 
54 K, respectively. The Cu-I boracite stays also cubic down 
to very low temperatures |7j and is expected to behave 
analogously to Cr-I boracite, since both Cr++ and Cu"'""'" 
are Jahn- Teller ions. However, the Neel temperature of 
Cu-I boracite is not known so far. Due to the cubic para- 
magnetic prototype point group 43ml', only two antifer- 
romagnetic point groups are allowed: 



4'3m' and 43m . 



(40) 



Whereas the point group 4'3m' allows the linear magneto- 
electric (ME) effect (invariant ctabEaHb) and both higher 
order ME effects (invariants ^^^^HaEbEc and P^^bcEaHtHc)., 
the point group 43m allows only the invariant P^^^EaH^Hc 
(see [52], Table 2). It could be that one or more of these 
three boracites have 4'3m' symmetry. 

In order to check this, ME measurements on single 
crystals would be desirable. However, the linear ME ef- 
fect can in principle also be detected on compacted, mag- 
netoelectrically annealed, polycrystalline powders or ce- 
ramics [^1231124] . both by MEh and MEe methods. The 



Pabc^aHbHc-eSect cancels out statistically in powders and 
ceramics because the sign of the coefficient is rigorously 
linked to the orientation of the non-centrosymmetric, nu- 
clear absolute structure (note, only in poled ferroelectric 
ceramics this statistics is broken and the effect becomes 
measurable in principle, when using the MEh effect). Be- 
cause the j^ij^HaEbEc-eSect cannot be measured by a 
MEh quasistatic method, it would not disturb when one 
measures the linear effect of a phase of the point group 
4'3m'. The j^bc^aEbEc-eSect can be measured by a MEe 
method or by a MEh technique in presence of a static 
electric field, as first shown by O'Dell [33]. As a conse- 
quence, even ME measurements on compacted powders 
or on ceramics would allow to distinguish between the 
two point groups, with 43m showing no ME effect at all 
and with 4'3m' only the linear effect when a MEh qua- 
sistatic method is used. We would like to suggest to grow 
single crystals [STllM] and to perform MEh measurements 
[3H] both on single crystals and polycrystalline samples. 
A polycrystalline probe of a cubic magnetoelectric crystal 
is obviously completely isotropic for the linear magneto- 
electric effect. Hence with a view on applications, the pro- 
duction of polycrystalline material could be simpler than 
growing single crystals. 



3.4.2 External magnetic octupole field of a cubic 
magnetoelectric crystal 

Dzyaloshinskii [13] hypothesized that there may exist an 
external magnetic field for free cubic (or axion type) mag- 
netoelectric media. Normally, the external magnetic field 
outside an antiferromagnet decays exponentially. Much 
earlier, Dzyaloshinskii [1^ predicted the existence of an 
external magnetic and electric field for magnetoelectric 
media with high symmetry, and specifically for the anti- 
ferromagnetic Cr2 03. Such a field decreases according to 
a power law from the surface of a body, in contrast to the 
exponential decay typical for substances with lower sym- 
metries. For Cr203, this external magnetic field behaves 
as the field of a magnetic quadrupole, which was subse- 
quently confirmed experimentally by Astrov et al. [5l[6]. 

For the magnetoelectric crystal with a cubic symme- 
try, which is characterized by the purely axion piece in 
the linear constitutive relation, we expect the similar ef- 
fect [13]. Taking, for definiteness, a ball made of such a 
material (recall that Astrov et al. in the 1960s used little 
balls made out of the antiferromagnetic Cr2 03), one can 
perform a usual expansion of the electric and magnetic 
fields into spherical harmonics. From a preliminary quali- 
tative analysis of the surface charges and currents induced, 
one can expect an external magnetic field of the octupole 
structure, namely. 



B 



(41) 



which is directly proportional to the axion piece a (here, 
as usual, the yim{0,'-p) denote the spherical harmonics). 
This preliminary result is based on the fact that is 
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invariant under the cubic symmetry group of the prospec- 
tive cubic axion materiaL 

3.5 Relations of a to gyrator, axion, and perfect 
electromagnetic conductor (PEMC) 

The constitutive law for the pseudoscalar piece, namely 
([5]) together with (fTS]) . reads 

Ha = -aEa , (42) 

see also (fT3| . (fT4| . This structure is not unprecedented, 
see the discussion in [20]. In electrical engineering, in the 
theory linear networks, more specifically in the theory of 
two ports (or four poles), Tellegen 64,65 invented the 
gyrator the defining relations of which are 

Vi ^ -3 12, 

V2= szi, (43) 

where v are voltages and i currents of the ports 1 and 2, 
respectively, see also [44] . Let us quote from Tellegen [6^ , 
p. 189: "The ideal gyrator has the property of 'gyrating' 
a current into a voltage, and vice versa. The coefficient 
s, which has the dimension of a resistance, we call the 
gyration resistance; 1/s we call the gyration conductance." 
The gyrator is a nonreciprocal network element that is 
discussed in the electrical engineering literature [ISllST] . 
for more recent developments see, e.g., [50] and [71] . 

For dimensional reasons, the electromagnetic excita- 
tions {D°',Ha) are related to the currents and the field 
strengths {Ea,B'^) to the voltages, that is, 

Ea = -sHa, 

= sD". (44) 

If we put s = 1/a, then ([44]l and ([42| coincide. Similar as 
the gyrator rotates currents into voltages, the axion piece 
'rotates' the excitations, modulo a resistance, into the field 
strengths. 

The next "isomorphic" structure to discuss is axion 
electrodynamics, see Ni (40j . Wilczek |70| . and, for more 
recent work, Itin [26ll27] l28] . We add to the usual Maxwell- 
Lorentz law for vacuum electrodynamics ©^"^ = Yq \/—gF^' 
an axion piece patterned after the last term in (jlOp . then 
we have the constitutive law for axion electrodynamics, 

^Yo^gF^'' + ]^ae'''"'''F^^. (45) 

In order to get a feeling for the Lagrangian of axion 
electrodynamics, we substitute ([T3]) . (I14p . and ([TS]) into 
the right-hand-side of ([4]); after some algebra we find 

^^^''Fxu ^^{B H D E) 

= l^B^ - lerSoE^ -a,Y„ E B . (46) 

^ V ^ axion part 

Maxwcll-Loiciitz part 



As we saw, in Cr203 we had ctr ~ 10"^. It is every- 
body's guess what it could be in elementary particle the- 
ory. There, in spite of extented searches via the coupling 
of photons to the hypothetical pseudoscalar axion field, 
amongst other methods, nothing has been found so far, 
see the review of Carosi et al. [TO] . 

These isomorphisms even found a further field of appli- 
cations: In 2005, Lindell & Sihvola [361137] . see also [35j, 
introduced the new concept of a perfect electromagnetic 
conductor (PEMC). It also obeys the constitutive law 
Some applications were studied by Jancewicz |30| and by 
Illahi and Naqvi [25] , see also the references given therein. 
The PEMC is a generalization of the perfect electric and 
the perfect magnetic conductor. In this sense, it is the 
'ideal' electromagnetic conductor that can be hopefully 
built by means of a suitable metamaterial, see Sihvola 
[57] . The pseudoscalar a is called Tellegen parameter by 
Lindell et al., see [38], p. 13 (for a more general view, see 
[58]): artificial Tellegen material has been produced and 
positively tested by Tretyakov et al. [S5], amongst others. 

Accordingly, (i) the gyrator, (ii) the axion field of ax- 
ion electrodynamics, (iii) the pseudoscalar a of Cr203, 
and (iv) the PEMC are isomorphic structures. Wilczek 
(private communication) remarked to these isomorphisms 
between these four systems: "It's a nice demonstration of 
the unity of physics." Already in his paper of 1987 [7^ he 
argued that "...it is... not beyond the realm of possibility 
that fields whose properties partially mimic those of ax- 
ion fields can be realized in condensed-matter systems." 
...This is, indeed, the case: The pseudoscalar a of Cr203 
is the structure Wilczek looked for. 

We would like to thank the organizers, in particular Nicola 
Spalding (Santa Barbara), for giving us the possibility to present 
a seminar at their workshop. We are grateful to Nicola Spalding 
for pointing out to us reference [5]. 
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